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We investigate experimentally and theoretically the be- 
havior of Coulomb blockade (CB) peaks in a magnetic field 
that couples principally to the ground-state spin (rather than 
the orbital moment) of a chaotic quantum dot. In the first 
part, we discuss numerically observed features in the magnetic 
field dependence of CB peak and spacings that unambiguously 
identify changes in spin S of each ground state for successive 
numbers of electrons on the dot, N. We next evaluate the 
probability that the ground state of the dot has a particu- 
lar spin S, as a function of the exchange strength, J, and 
external magnetic field, B. In the second part, we describe 
recent experiments on gate-defined GaAs quantum dots in 
which Coulomb peak motion and spacing are measured as a 
function of in-plane magnetic field, allowing changes in spin 
between TV and N + 1 electron ground states to be inferred. 

I. INTRODUCTION 

In the absence of interactions, electrons populate the 
orbital states of a quantum dot or metallic grain in an 
alternating sequence of spin-up and spin-down states, in 
accordance with the Pauli picture. In this case, the total 
spin of the dot is either zero, when the number of elec- 
trons N is even, or one half at odd N. No higher spin 
states appear. It is well known that in the presence of 
interactions, this simple scheme can be violated. Such 
deviations from simple even-odd filling may appear as 
Hund's rules in atomic physics, or be as extreme as com- 
plete spin polarization leading to ferromagnetism by the 
mechanism of the Stoner instability. 

By measuring electron transport through a weakly cou- 
pled quantum dot at low temperature and bias voltage 
(compared to the quantum level spacing) one can effec- 
tively study differences in ground state (GS) properties of 
the dot at successive electron number, providing a means 
of investigating the actual filling scheme Q. A stan- 
dard experimental approach Q is to measure conduction 
through the dot via two tunneling leads as a function of 
a voltage, V g , applied between the dot and a third elec- 
trically isolated electrode, known as a gate. Because the 
dot GS energy, En, at a fixed electron number is influ- 
enced by V g , the gate can be used to set the number of 
electrons on the dot. The fact that a large Coulomb en- 
ergy is needed to add a single electron to the dot typically 
suppresses conduction through the dot in the tunneling 



regime, where whole charges must tunnel for conduction 
to occur. This effect is known as the Coulomb blockade 
(CB). However, at specific values of gate voltage, denoted 
Vg N \ where the condition 25jv_i(Vg) = E^(V g ) is satis- 
fied, conductance increases dramatically, clearly marking 
this degeneracy condition in an experimental trace. The 
position, Vj, , of the N th peak in the conductance is pro- 
portional to //at = En(0) — En-i(0). Accordingly, the 
distance between successive peaks is = /zjy — /xjv-i- 
(We take the constant of proportionality converting gate 
voltage to dot energy to be unity for the theoretical dis- 
cussion; experimentally, this constant can be readily mea- 
sured, for instance, by comparing the influence of bias 
and gate voltages.) 

In the last few years the GS spin of a variety of nanos- 
tructures, including metallic grains Q, semiconducting 
quantum dots ||-|(|, and carbon nanotubes j7|||, have 
been investigated using CB peak motion in a magnetic 
field. If one neglects the magnetic field coupling to the 
orbital degrees of freedom, one expects the field to mani- 
fest itself only through the Zeeman splitting, resulting 
in a shift of the GS energy by gfisSB, where fiB is 
the Bohr magneton and S denotes the GS total spin. 
For 2D quantum dots, orbital coupling can be strongly 
suppressed — though not eliminated entirely — by orient- 
ing the field strictly in the 2D plane. This is the exper- 
imental approach that will be described in Section III 
of this paper. On the other hand, for ultrasmall grains 
and nanotubes, it is reasonable to ignore orbital coupling 
for any field direction since for practical magnetic fields 
the total flux cutting the structure is much less than the 
quantum of flux. 

For semiconducting quantum dots ||-|(| the Zeeman 
splitting for S — 1/2 becomes comparable to the mean 
single electron level spacing 

Si = (s*); Si = Ei - e,+i, (1) 

(here £j denotes the orbital energy of the one-electron 
orbital state i, and (...) stands for an average over dif- 
ferent levels) at B ~ 1 T for the GaAs dots is Refs. [g||] 
and B ~ 10 T for the Si dot in Ref. g|. At lower fields, 
the splitting seldom exceeds Sj, and one might expect a 
simple even/odd Pauli picture. In this case, CB peak po- 
sitions would move with B as consecutive pairs, creating 
a pattern of alternating downward and upward moving 
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peaks, as seen in Fig. |l|a. The trajectories of the pair of 
peaks would form straight lines with the same slope of 
g/ib/2- The peak spacings, A^, will have a similar pat- 
tern with slopes of g/ib- Spin-orbit interaction can lead 
to fluctuations in the g factor, resulting in fluctuations 
in the slopes of the lines, J9|,[l0|], but will not change the 
pattern of downward and upward moving peaks. Once 
g/ibB exceeds a particular spacing Sj, the peaks should 
cross. 

Experimentally, the Pauli picture appears to work well 
for metallic grains while for carbon nanotubes the 
situation is less clear: one published experiment finds 
simple even/odd filling [Q, and one finds a more com- 
plicated scheme §]. For multielectron semiconducting 
dots HUD the data are somewhat confusing, but ap- 
pear to disagree with the simple even/odd filling scheme 
||||]|]ll|] ' presumably as a result of electron-electron in- 
teractions. In symmetric, few-electron vertical dots 
one may again recover relatively simple behavior for the 
first few electrons, with a well-understood appearance 
of Hund's rules. These conclusions are based on exper- 
iments where orbital magnetism dominates spin. How- 
ever, for few electrons (i.e., TV < 10) the filling scheme 
can be readily interpreted, allowing the full spin structure 
to be mapped out as a function of N and B @ . 

Figure |l| shows CB peak motion and spacing as a func- 
tion of magnetic field (coupling only to spin) for both a 
nonintcracting electron model and an interacting elec- 
tron model presented in detail in Section II. The behav- 
ior of this model appears more similar to the experimen- 
tally observed behavior in semiconducting quantum dots 
(see Fig. ^|). In some cases peak positions (/ijv) move as 
function of the magnetic field in bunches of two in the 
same direction and with the same slope, while their cor- 
responding spacings (A^) do not change as function of 
the magnetic field. Moreover, not every change in the 
direction of the motion of a peak can be explained as a 
crossing of two orbital levels. Thus, a description of the 
system beyond the Pauli picture is needed. 

We will present a theory that goes beyond the simple 
Pauli picture and explains the puzzling features of the 
magnetic field dependence of the conductance peaks po- 
sitions as a result of spontaneous spin polarization of the 
electrons in the dot ju). The possibility of spontaneous 
magnetization in quantum dots has been considered in 
different contexts [[iTj-pof . For instance, it has also been 
suggested as an explanation |li| for the absence of a bi- 
modal distribution in the conductance peak spacings pl| . 
Such a magnetization was also predicted to lead to kinks 
in the parametric motion of the peaks, (e.g., due to or- 
bital effects of a perpendicular magnetic field) Mj. 
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FIG. 1. Typical conductance dependence on the gate volt- 
age and magnetic field, B, with and without exchange in- 
teraction: (a) J = 0, (b) J = 0. l<5i . [i' N is determined as 
Mjv = MJV — e 2 /C and taken in arbitrary units. for each 
N is shifted arbitrarly by Nx constant to avoid overlapping. 



II. THEORY 

We begin the theoretical discussion by reviewing 
the considerations determining the general form of the 
Hamiltonian describing the properties of a chaotic dot in 
the metallic regime. A more complete description may 
be found in Rcf. |Q. The single electron spectrum 
is characterized by the mean level spacing 5i and the 
Thouless energy Et ~ h/t erg , where t erg is the time it 
takes for a classical electron to cover the energy shell in 
the single-particle phase space. For diffusive and ballis- 
tic systems t erg equals to L 2 /D and L/vp respectively, 
where Vf is the electronic Fermi velocity, L is the system 
size, and D is the diffusion coefficient. The dimension- 
less conductance, gx — Ex /Si, in the metallic regime 
is large, i.e., jf > 1. In this regime the statistics of 
the single electron spectrum on scales smaller than Et 
are well described by random matrix theory (RMT) [§2| , 
which gives a quantitative description of the phenomenon 
of level repulsion. For an ensemble of M x M matrices 
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(M — > oo) with random and independent elements, the 
probability density of a realization of the spectrum Si is 
given by [pi: 



P{e,i) oc exp 



2 £f 



5l 



(2) 



where /3 is equal to 1, 2 or 4 for the orthogonal, unitary 
and symplectic ensembles respectively. The orthogonal 
(unitary) RM ensemble corresponds to weakly disordered 
dots with preserved (violated) time reversal symmetry 
and negligible spin-orbit. We assume that the dot is small 
enough for the spin-orbit interaction to be negligible and 
thus avoid discussion of the symplectic ensemble (/3 = 4) 
& 

The second step is to consider effects of electron- 
electron interactions. For the simplest case of short-range 
interactions, we write 



H int {r) = XtxV&lf), 



(3) 



where V oc L d is the volume of the dot (in d dimensions) 
and A is a dimensionless coupling constant characteriz- 
ing the strength of the interaction. Matrix elements of 
this interaction, in the basis of eigenstates ipi(r) of the 
noninteracting Hamiltonian, are given by 



M% = \8xV 



drip* (r)ip* (r)ip k (r)ipi (r) . 



(4) 



It is important to note that the statistical properties of 
the interaction matrix elements are completely de- 
termined by the statistical properties of the single elec- 
tron eigenstates ipi(r) and cannot be chosen in an arbi- 
trary fashion as in the random interaction model p0[ . 
In fact, in the limit gx — ► oo the diagonal matrix el- 
ements (i,j,k,l pairwisc equal) are self-averaged by the 
space integration, Eq. ([|) , and thus show no level-to-level 
or sample-to-sample fluctuations |l4|. In the same limit 
the off-diagonal matrix elements Eq. (^ turn out to be 
negligable. 

Using the statistical properties of the interaction ma- 
trix elements it turns out that a large class of disordered 
metallic dots can, under very general conditions, be de- 
scribed by a remarkably simple Hamiltonian with only 
three coupling constants, which do not fluctuate. In the 
limit of large gx, the interaction part of the Hamiltonian 
corresponds to 



H mt = E C N 2 - J{S f + X BCS T ] f, 



(5) 



(terms linear in N are allowed, but they can be included 
into the one-particle part of the Hamiltonian) where 
N is the number operator, S is the spin operator and 
T = X)i c iT c *l ( C *T annihilates an electron in the i th sin- 
gle electron orbital with spin f). In the simple model 
with short-range interaction and preserved time reversal 



((3 = 1) the above coupling constants have the following 
form: 

E c = ^XSx; J = 2XSr, Xbcs = XS V (6) 

If time invariance is broken (/? = 2), Xbcs = since the 
operator T is incompatible with the symmetry. Note that 
the interaction Hamiltonian Eq. ([}]) represents a particu- 
lar model, and the expressions for the coupling constants 
Eq. (^) are valid only for this model. At the same time, 
the effective interaction Hamiltonian, Eq. (|J) , is more ro- 
bust and depends only on the symmetries of the problem 
(for instance, on the absence of the spin-orbit scattering) 
and on the condition gx 3> 1. 

The first two terms in Eq. represent the dependence 
of the energy of the dot on the total number of the elec- 
trons and on the total spin respectively. They commute 
with each other and with the single-particle part of the 
Hamiltonian. Therefore, all states of the grain can be 
classified by N and S. The term proportional to Xbcs 
appears only in the orthogonal case ((3 = 1). Provided 
that Xbcs < this term leads to the superconducting 
instability. Superconducting correlations are suppressed 
by the magnetic field, and thus do not exist for /3 = 2. 

Thus, the general form of the Hamiltonian describing 
electrons in a chaotic dot which is not superconducting 
is given by 

H = £ * n * + E C N 2 - JS(S + 1) + gUbSB. (7) 

i 

Note that the only random component of the problem 
is the single-particle spectrum et, while the exchange J 
and charging energy E c do not fluctuate. It should be 
realized, though, that Eqs. (||J^) become exact only in 
the limit gx — ► oo. The corrections to Eqs. (HQ), which 
appear at finite gx are sometimes important. However, 
if gx 3> 1, these corrections do not bring essentially new 
physics to the problem of spin magnetization. 

The conductance is calculated using the many-particle 
energies En and wave functions obtained numerically for 
GOE and GUE random matrix realizations. An example 
of the peak positions and peak spacing evolution as func- 
tions of the magnetic field for a particular GOE realiza- 
tion is presented in Fig. [l]for (a) J = and (b) J = 0.15i. 
For both cases we present [x' N — //jy — E c . The noninter- 
acting case, Fig. [l|a, shows all the previously described 
features of the Pauli behavior. Once a weak exchange in- 
teraction is included the behavior changes qualitatively, 
not just at high field, but also in the vicinity of B = 0. In 
particular, peak positions are not always paired, result- 
ing in occasional field- independent peak spacings. This 
occurs when two consecutive orbitals are first filled with 
down spin electrons and only later they acquire up elec- 
trons. Generally the enhancement of the spin of the dot 
by S will be accompanied by a 2S bunch peaks moving 
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with the same slope. If the peak spacing is plotted, 
two sets of 25 — 1 flat curves sandwiching a sloped one 
will appear. Changes to the crossing also attributed to 
a sudden change of the GS spin which is not associated 
with a single electron orbital crossing. 

As illustrated above, the appearance of higher spin 
clearly manifests itself in the peak position and peak 
spacing trajectories. Using Eq.(^) one can predict the fre- 
quency of spontaneous magnetization appearances. For 
weak exchange, J <C 5% , the probability that the dot GS 
has a spin 5 is determined by the probability of finding 
25 orbitals so close to each other that the gain in ex- 
change energy (due to the polarization) overwhelms the 
loss in the kinetic energy (due to single rather than dou- 
ble occupations of the orbitals) . 

Thus the probability for a certain value of ground state 
magnetization 5 boils down to the question of the prob- 
ability of finding 25 close single electron levels. That 
probability may be estimated using random matrix the- 
ory, as we now present. 

The probability of finding a set of single electron or- 
bital at energies £, is given by Eq. (^). Thus as long 
as 



£%+2S - £» < JS(S + 1) + gusBS, 



(8) 



a ground state spin of at least 5 will appear. Taking into 
account the probability of Eq. (||), it is possible to write 
the probability of obtaining a ground state spin 5 as a 
function of the linear combination: 



9VbB 



(9) 



Si IS + I\5, 
resulting in the following probability 

P^ B (S) = c|jf 0W+i)(2S-U(i - K^X 2 ), (10) 

The coefficients Cg and Kg depend on both /3 and 5. 
Their values for 5 = 1,3/2 are presented in table I. 
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TABLE I. The factors C§ and K§ appearing in Eq. ((Toj) . 

Thus the RMT model presents a striking conclusion. 
The influence of exchange interactions and externally ap- 
plied magnetic field on the ground state magnetization of 
a quantum dot may be summed up in a single parameter 
scaling function, where the scaling parameter is a simple 
linear combination of J and B. As has been shown from 
numerical simulations |23|, the scaling holds for larger 



One may wonder how well does the above picture hold 
for systems for which real correlation between the elec- 
tron exist, and the dimensionless conductance is not too 
large. A canonical example for such systems is the Hub- 
bard model: 



+ U X! a !,T ai >T a !,i a U + 9VbSB (11) 



where < i,j > denotes nearest neighbor lattice site, a\ a 
is an creation operator of an electron at site i with spin <r, 
ti is the site energy, chosen randomly between —W/2 and 
W/2 with uniform probability, and U is the interaction 
constant. 

The Hilbert space is large even for relatively small sys- 
tems. For example for a 4 x 4 lattice with 6 electrons 
for the 5 = sector has the size of 313600. Using the 
Lanczos method we obtain the many-particle eigenval- 
ues En(S) as function of the ground state in each spin 
sector 5 for different values of interaction U . The dis- 
order was chosen as W = 8t, which corresponds to the 
metallic regime, although the value of the dimensionless 
conductance is quite low qt ~ 0(1). 

We calculated the probability for the appearance of a 
specific value of GS spin 5 for different values of U and 
magnetic field B. In order to use the scaling parameter 
X we need to deduce the appropriate value of J for each 
value of U. This has been done by fitting the dependence 
of the average lowest energy in a given spin sector (E(S)) 
to 



(E(S)) =S lS 2 - J5(5+l). 



(12) 



values of X than expected from the first and second term 
perturbative analysis presented in Eq. (JIG). 



The distribution Pj(u),b{S) based on numerical results 
for 1000 different realizations of the Hubbard model, for 
different values of U and B for /3 = 1 and f3 = 2 is pre- 
sented in Fig. H Although some deviations do appear, 
especially for the higher values of U, one can see that 
the overall form of the scaling function holds remarkably 
well. Thus, though the Hamiltonian in Eq. (Q) is obtained 
under the conditions of large qt and no correlations in 
the system, it nonetheless appears to capture much of 
the physics even for a moderate gx and in the presence 
of electron-electron correlations. The reason for this is 
that although neither J nor E c is constant once qt ~ 1, 
their fluctuations remain small relative to their average, 
as is demonstrated in Fig. ||. 



III. EXPERIMENT 

In this section we describe recent experimental mea- 
surements of GS spin for a gate-defined GaAs quantum 
dot containing roughly 400 electrons. As discussed in the 
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FIG. 2. Probabilities of S — 1 for J = 0.1<5i as functions of 
the magnetic field B for (a) the GOE and (b) GUE cases. The 
small symbols represent numerical results for the Hubbard 
model, the large filled symbols represent numerical results for 
the RMT model, while the curves (solid - first order, dashed 
- second order) represent Eq. (0). 

Introduction, the dot is coupled to electron reservoirs via 
tunnelling leads (i.e., g < 2e 2 /h for both leads) so trans- 
port is dominated by CB effects. Measurements were 
carried out at sufficiently low temperature and bias that 
the differences in GS energies between dots with N + 1 
and N electrons should be extractable from CB peak po- 
sitions, Vg. 

The dot is formed at the interface of a GaAs/AlGaAs 
heterostructure (90 nm below the wafer surface) by elec- 
trostatic depletion using surface gates. A Si delta dop- 
ing region is located 40 nm above the heterointerface. 
The two-dimensional electron gas (2DEG) has density 



2.0 xlO 11 . 



and bulk mobility ~ 1.4 X 10 cm /Vs, 



yielding a transport mean free path ~ 1.5 nm. The 
small dot area, A ~ 0.25 /im 2 , makes transport predom- 
inantly ballistic within the device. Characteristic en- 
ergy scales for the measured device include the mean 
level spacing, A = 2irh 2 /m*A ~ 30 fieV, the charg- 
ing energy, Eq ~ 400ueV, and the Thouless energy 
E t h = TivpA^ 1 / 2 ~ 340/ieV. Measurements were carried 
out in a dilution refrigerator with a mixing chamber tem- 
perature of 25 mK using standard ac lock-in techniques 
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FIG. 3. The exchange parameter J as functions U for a 
4x4 Hubbard model with 6 electrons. 



with a source-drain bias voltage of 2 ueV. A base elec- 
tron temperature of T e ~ 50 mK was determined from 
CB peak widths. 

To allow the magnetic field to couple predominantly to 
spin, the sample was oriented with the plane of the elec- 
tron gas along the axis of the primary solenoid, aligned 
manually to within 0.5 degrees. In addition, a small pair 
of coils attached to the vacuum can of the fridge, ori- 
ented perpendicular to the plane of the sample (see Fig. 
4), was used to null out any perpendicular field from 
misalignment as well as to explicitly break time-reversal 
symmetry. Both primary and trimming coils were un- 
der computer control, allowing sweeps of strictly paral- 
lel field by trimming out any perpendicular component 
for any applied parallel field. (We estimate the uncer- 
tainty in B±_ to be less than o /4 through the dot at 
B\\ = 5T.) Despite this precise field trimming capabil- 
ity, similar measurements in larger dots, fabricated 
on the same wafer, indicate orbital coupling due to the 
strictly parallel field for B\\ > 0.5T (evident, for instance, 
from the disappearance of the weak localization feature 
at B±_ — at higher Bn). The origin of this surprising 
coupling — which makes interpretation of our high-field 
data difficult — is still under investigation. 

Conductance measurements across ten consecutive 
Coulomb blockade peaks, measured as a function of V g 
and Sii (i.e., strictly Bu , properly trimmed), are shown in 
Fig. [|a. More positive gate voltage corresponds to higher 
energy, and can be calibrated from the CB "diamonds" 
using high source-drain bias measurements pj. All data 
are taken with B±_ = 2§mT in order to ensure that time- 
reversal symmetry is broken, which changes the statistics 
of dot wave functions. 

In addition to the individual motions of each peak, 
there is a diamagnetic shift of all peaks (see Fig. ||b), 
presumably due to the effect of the parallel field on the 
effective well confinement potential. The slight param- 
agnetic shift visible at low field, (Bh < 0.2T), is not 
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FIG. 4. a) Schematic diagram of a gate-defined quantum 
dot, indicating the orientations of B11 and B± relative to the 
planar quantum dot. (Orientation of Bu within the plane is 
not accurately depicted and has not been investigated.) b) 
Diagram showing placement of superconducting coils used to 
generate B± on vacuum can, inside primary solenoid used to 
generate Bm . 

understood at present. In the following analysis of data, 
the common curve (in Fig. ||b) has been subtracted from 
each peak position. 

Peak positions and peak spacings extracted from the 
peaks in Fig. 5 are shown in Fig. 6. The slope of peak 
positions as a function of Bu is consistent with a Zeeman 
energy term Es = ±-^g/.iBB, using the g-factor for bulk 
GaAs, \g\ — 0.44. As discussed in the Introduction, alter- 
nating slopes for consecutive peaks would indicate an al- 
ternating 0, i, 0, i... GS spin structure. Our data, on the 
other hand, shows three consecutive pairs of peaks mov- 
ing with the same slope, suggesting the presence of higher 
spin states. Proposed values for the eight consecutive GS 
spin states shown here are included in Fig. 6a. We em- 
phasize, however, that these are only plausible values for 
the spin; it is not possible to determine unambiguously 
the absolute magnitude of GS spin from measurements 
of peak position, which reflect changes in spin from the 
N to N+l ground states. The values shown in Fig. 6a 
minimize the ground state spin for the system. 

In the proposed spin labelling scheme, three out of the 
five even-iV states have 5=1, i.e., P(S = 1) ~ 0.6. 
This fraction of S = 1 relative to the number of S = 
states is well beyond the expected value given reasonable 
estimates of J for GaAs dots. Of course, with only five 
spin states considered, statistics are quite poor. Further 
experiments are needed to see if this discrepancy is sig- 
nificant. 

Figure 7 shows that peak spacings clearly separate into 
three branches, a top branch with slope roughly gfi (cor- 
responding to a GS spin decrement followed by an in- 
crement) a bottom branch with slope roughly — g\i (cor- 
responding to a GS spin increment followed by a decre- 
ment) and a middle branch with slope near zero (corre- 
sponding to two consecutive increments or decrements). 
The existence of the middle branch is the signature of 
higher GS spins. The good agreement between the slopes 
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FIG. 5. a) Nine consecutive Coulomb blockade peaks mea- 
sured as a function of gate voltage V g and parallel field B. 
Conductance is shown in grayscale, with black indicating high 
conductance and white low conductance. The field shown on 
the bottom axis is the field strictly parallel to the heterointer- 
face; perpendicular field is held constant at 20mT as discussed 
in the text. Note that conductance in the valleys is low but 
not zero, due to strong tunnelling in the leads, b) Average 
peak position of the peaks in a), showing a diamagnetic shift 
to higher energies common to all peaks. This average is sub- 
tracted from all peak positions before further analysis. 

of the upper and lower branches and the expected slopes 
of g/j>, as well as the absence of a range of intermedi- 
ate slopes suggest that the peak spacing reflects spin 
rather than orbital coupling. At higher fields, the di- 
rections of peak motion change, often abruptly and from 
one straight segment to another, as seen in Fig. 8. This 
behavior is qualitatively similar to the numerical data in 
Fig. lb. The rounding of straight segments where the 
slope changes presumably results from spin-orbit inter- 
action which mixes spins, and may provide a direct mea- 
sure of spin-orbit interactions in dots. It is interesting to 
note that the peaks analyzed in Figs. 6 through 8 were 
measured in a regime of high tunneling conduction in the 
leads. This can be seen by noting the grayscale of Fig. 5. 
When the dot is more pinched off from the reservoirs, so 
that the CB peaks have a height of 0.1 e 2 /h or less, peak 
motion is more difficult to interpret, and does not seem 
to follow the clear patterns illustrated, for instance, in 
Fig. 8. As seen in Fig. 9, peak spacings for four consecu- 
tive weak-tunneling peaks at do not move with slopes of 
5/i at low Bu. Though the data are noisier in this con- 
figuration, it is clear that the slopes are in all cases less 
than g/i. We do not have an explanation for the different 
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FIG. 6. a) Peak positions (with average motion subtracted) 
as a function of parallel field for the nine peaks shown in Fig. 
5a. Straight lines indicate expected peak motion for spin-i 
transitions using the g factor for bulk GaAs, \g\ = 0.44. Pos- 
itive slopes indicate transitions to lower spin states; negative 
slopes indicate transitions to higher spin states. Numbers 
indicate a possible ground-state spin structure for ten consec- 
utive ground states. Peak position data are offset for clarity, 
b) Peak spacings for the position data shown in a), offset for 
clarity. Solid lines indicate expected motion of peak spacing 
for spin- g transitions, using \g\ = 0.44 



behavior depending on if the dot is more closed off or 
nearly open. 



IV. DISCUSSION 

In conclusion, we have shown that a relatively weak ex- 
change interaction qualitatively explains the deviations 
from the Pauli picture seen in recent experiments ||f|,D ■ 
We calculate the probability that different values of spin 
appear as a function of exchange interaction and mag- 
netic field coupling to the spin. The strength of the 
exchange interaction, J, and the effective g factor are 
the only adjustable parameters that determine the prob- 
abilities for a GS of the dot to have any given spin at 
any given magnetic field for both GOE and GUE cases. 
In particular we predict that these probabilities will fol- 
low a one-parameter scaling law over a wide range of 
magnetic fields and exchange interaction strengths. We 




FIG. 7. Peak spacings from Fig. 6b, offset to align spacings 
at By = 0. Slopes cluster in three branches, interpreted as a 
change from decreasing to increasing spin transitions (positive 
slope), two consecutive spin transitions in the same direction 
(zero slope), and a change from increasing to decreasing spin 
transitions (negative slope). Solid lines indicate expected mo- 
tion of peak spacing for spin-i transitions, using \g\ = 0.44, 
with no adjustable parameters. 



also report preliminary experiments to investigate the ex- 
change effects in the ground state spin of quantum dots. 
The signature of interaction effects is the appearance of 
higher-spin ground states, which show up in the experi- 
mental data as peak spacing traces that have zero slope 
as a function of parallel magnetic field. Such features are 
indeed observed, as seen for instance in Fig. 8. Further 
experiments are needed to obtain sufficient statistics to 
quantitatively test the predictions of theory. 
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